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Abstract 



The properties of high-mass multijet events produced at the Fermilab proton- antiproton 
coUider are compared with leading order QCD matrix element predictions, QCD parton 
shower Monte Carlo predictions, and the predictions from a model in which events are 
distributed uniformly over the available multibody phase-space. Multijet distributions 
corresponding to (4N-4) variables that span the N-body parameter space are found to 
be well described by the QCD calculations for inclusive three-jet, four-jet, and five-jet 
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events. The agreement between data, QCD Matrix Element calculations, and QCD 
parton shower Monte Carlo predictions suggests that 2 — > 2 scattering plus gluon ra- 
diation provides a good first approximation to the full LO QCD matrix element for 
events with three, four, or even five jets in the final state. 
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1 Introduction 



A study of the properties of events containing three-or-more jets produced in high- 
energy hadron-hadron coUisions can (i) test our understanding of the higher-order QCD 
processes that result in multijet production, (ii) test the QCD parton shower Monte 
Carlo approximation to the full leading order (LO) QCD matrix element, and (iii) 
enable a search for new phenomena associated with the presence of many hard partons 
in the final state. The CDF collaboration has recently reported the characteristics 
of two-jet, three-jet, four-jet, five-jet, and six-jet events ^ produced at the Tevatron 
proton-antiproton collider operating at a center-of-mass energy of 1.8 TeV. Results from 
an analysis of events with multijet masses exceeding 600 GeV/c^ were presented for a 
data sample corresponding to an integrated luminosity of 35 pb^^. The multijet-mass 
distributions, leading-jet angular distributions, and mass dependent jet multiplicity 
distributions were shown to be well described by both the NJETS 0] LO QCD matrix 
element calculation for events with up to five jets, and the HERWIG P| QCD parton 
shower Monte Carlo calculation for events with up to six jets. For these selected 
distributions the QCD predictions were found to give a good description of the data. 

In the present paper we use a larger data sample and a more comprehensive set of 
multijet distributions to extend our comparison of the properties of high-mass multijet 
events with QCD predictions. In particular, we use the set of (4N-4) variables that 
span the N-jet parameter space and were recently proposed by Geer and Asakawa 0], 
and compare the observed three-jet, four-jet, and five-jet event characteristics with (a) 
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NJETS LO QCD matrix element predictions, (b) HERWIG parton shower Monte Carlo 
predictions, and (c) predictions from a model in which events are uniformly distributed 
over the available multijet phase-space. Results are based on a data sample which was 
recorded by the CDF collaboration during the period 1992 - 1995, and corresponds to 
an integrated luminosity of 105 pb~^. 

2 Experimental Details 

A description of the CDF detector can be found in Ref. [^. Full details of the CDF 
jet algorithm, jet corrections, and jet resolution functions can be found in Ref. 0, 
and a description of the trigger and event selection requirements for the high-mass 
multijet sample are given in Ref. [0]. In the following we give a summary of the main 
details of the CDF detector, jet reconstruction, and event selection requirements that 
are relevant to results presented in this paper. We use the CDF co-ordinate system in 
which the origin is at the center of the detector, the z-axis is along the beam direction, 
9 is the polar angle with respect to the z-axis, and is the azimuthal angle measured 
around the beam direction. 

The multijet analysis described in the following sections exploits the CDF calorime- 
ters, which cover the pseudorapidity region |?7| < 4.2, where rj = — ln(tan^/2). The 
calorimeters are constructed in a tower geometry in r/ - space, and are segmented in 
depth into electromagnetic and hadronic compartments. The calorimeter towers are 
0.1 units wide in t]. The tower widths in are 15° in the central region and 5° at larger 
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|?7| (approximately |?7| > 1.2). Jets are reconstructed using an algorithm that forms 
clusters from localized energy depositions in the calorimeter towers. Calorimeter tow- 
ers are associated with a jet if their separation from the jet axis in [r], 0)-space AR — 
(A?7^+ A0^)^/^ < Rq. For the analysis described in this paper the clustering cone radius 
it!o = 0.7 was chosen. With this Rq a plot of the separation between all jets observed in 
the data sample described below reveals that, to a good approximation, clusters with 
separations AR < 0.8 are always merged by the jet algorithm into a single jet, and 
clusters with separations AR > 1.0 are never merged. Thus, the effective minimum 
observable separation between jets ARmin — 0.9 ± 0.1. Jet energies are corrected 
for calorimeter non-linearities, energy lost in uninstrumented regions and outside of 
the clustering cone, and energy gained from the underlying event. The jet corrections 
typically increase jet energies by 25% for jets with transverse energy Et — E sin 9 > 
60 GeV, where 9 is the angle between the jet axis and the beam direction. The jet 
corrections are larger for lower Et jets, and typically increase jet energies by about 
30% (40%) for jets with Et = 40 GeV (20 GeV). After correction, jet energies are 
measured with a precision as/ E of approximately 0.1 and multijet masses calculated 
from the jet four- vectors are measured with a precision (Tm/m of approximately 0.1. 
The systematic uncertainty on the jet energy scale is 5% for jets in the central region 
(|?7| < 1.2). There is an additional systematic uncertainty of 2% on the energy scale of 
jets at larger \ri\ relative to the corresponding scale for central jets. 

The data were recorded using a trigger which required Et > 300 GeV, where the 



9 



sum is over the transverse energies (Et) of all uncorrected jets with Et > 10 GeV, 
and the jet transverse energies were calculated assuming an event vertex at the center 
of the detector (x=y=z=0). In the subsequent analysis the J2Et was recalculated 
using the reconstructed vertex position and corrected jet energies, and summing over 
all jets with corrected Et > 20 GeV. Events were retained with Et > 420 GeV. To 
reject backgrounds from cosmic ray interactions, beam halo, and detector malfunctions, 
the events were required to have (i) total energy less than 2000 GeV, (ii) a primary 
vertex reconstructed with \z\ < 60 cm, (iii) no significant energy deposited in the 
hadron calorimeters out-of-time with the proton-antiproton collision, and (iv) missing- 
Et significance ^ = /(E^t Y^"^ < 6, where ^t = \ EE^il, and E^i is a vector 
that points from the interaction vertex to calorimeter cell i and has magnitude equal 
to the cell Et- These requirements select 30245 events. 



3 QCD and Phase-Space Predictions 

In the following we will compare observed multijet distributions with predictions from 
(a) the HERWIG g QCD parton shower Monte Carlo program, (b) the NJETS [| LO 
QCD 2 N matrix element Monte Carlo program, and (c) a model in which events 
are distributed uniformly over the available N-body phase-space. 
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3.1 The HERWIG Parton Shower Monte Carlo Calculation 

HERWIG 0] is a QCD parton shower Monte Carlo program that includes both initial- 
and final-state gluon radiation. HERWIG predictions can be thought of as QCD 2^2 
predictions with gluon radiation, color coherence, hadronization, and an underlying 
event. We have used version 5.6 of the HERWIG Monte Carlo program together with 
a full simulation of the CDF detector response. In our HERWIG calculations we 
have used the CTEQIM structure functions and the scale = stu/2(s^+u^+t^). 
HERWIG generates 2^2 processes above a specified p^"^*^ where p^"^*^ is the pt of 
the outgoing partons from the hard scatter before any radiation has occurred. It is 
important to chose a low value of p^"'"'^ so that adequate account is taken of events 
in which the detector response has fiuctuated upwards by several standard deviations 
and/or the spectator system accompanying the hard scattering process, including the 
initial state radiation, makes an unusually large contribution to the Et ■ We have 
set the minimum p^"^'^ to 60 GeV/c. The contribution to the J2Et > 420 GeV Monte 
Carlo sample from events with pl^^'^ < 60 GeV/c is estimated to be less than 1%. 

3.2 The NJETS QCD Matrix Element Calculation 

The NJETS Monte Carlo program provides parton-level predictions based on the 
LO QCD 2 —>■ N matrix elements. The calculation requires the minimum separa- 
tion between the final state partons in (77, 0)-space to exceed ARmin- We have set 
^Rmin = 0.9, and have used the KMRSDO structure function parameterization 
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with the renormahzation scale chosen to be the average pr of the outgoing partons. 
NJETS does not use a parton fragmentation model. Jet definitions and selection cuts 
are therefore applied to the final state partons. To enable a direct comparison between 
NJETS predictions and observed distributions we have smeared the final state parton 
energies in our NJETS calculations with the Gaussian jet energy resolution function: 

aE = 0.1 E . (1) 

This provides a good approximation to the CDF jet resolution function. 

3.3 Phase-Space Model 

We have used the GENBOD phase-space generator [|l^] to generate samples of Monte 
Carlo events for which the multijet systems uniformly populate the N-body phase- 
space. These phase-space Monte Carlo events were generated with single-jet masses 
distributed according to the single-jet mass distribution predicted by the HERWIG 
Monte Carlo program. In addition, the multijet mass distributions were generated 
according to the corresponding distributions obtained from the HERWIG Monte Carlo 
calculation. Comparisons between the resulting phase-space model distributions and 
the corresponding HERWIG and NJETS Monte Carlo distributions help us to under- 
stand which multijet parameters are most sensitive to the behaviour of QCD multijet 
matrix elements. 
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4 Multijet Variables 

To completely specify a system of N jets in the N-jet rest-frame we require (4N-3) 
independent parameters. However, the N-jet system can be rotated about the beam 
direction without losing any interesting information. Hence we need only specify (4N- 
4) parameters. We will use the N-jet mass and the (4N-5) dimensionless variables 
introduced and discussed in Ref . . In the following the variables are briefly reviewed. 

4.1 Three- Jet Variables 

In previous three-jet analyses ^, |12| it has become traditional to label the incoming 



interacting partons 1 and 2, and the outgoing jets 3, 4, and 5, with the jets ordered 
such that E3 > E4 > E^, where Ej is the energy of jet j in the three-body rest-frame. 
At fixed three-jet mass rrisj a system of three massless jets can be specified in the 
three-jet rest-frame using four dimensionless variables, X3, X4, cos6'3, and ip^, which 
are defined: 

(i) The Dalitz variables X3 and X4: 



2 Ei , , 



(ii) The cosine of the leading jet scattering angle: 



'Pav ■ 



cos ^3 = , (3) 

I Pav 1 1 -t 3 1 
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where the average beam direction: 

Av = (4) 

and particle 1 is the incoming interacting parton with the highest energy in the 
laboratory frame. 



(iii) '03, defined in the three-jet rest-frame as the angle between the three-jet plane 
and the plane containing the leading-jet and the average beam direction: 



_ P3X Pav)-{P,^ P5 ... 
P3 X PavWPa X F5 



To specify a system of three massive jets we must supplement the traditional three-jet 
variables with three additional parameters that describe the single-jet masses. These 
parameters are taken to be the single-jet mass fractions /a, /4, and /s, where: 



Thus we have eight three-jet variables which consist of maj, four parameters that 
specify the three-jet configuration (X3, X4, cos^a, and ■^a), and three variables that 
specify the single-jet masses {fa, f^, and /g). 
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4.2 Four- Jet and Five- Jet Variables 

A multijet system with more than three-jets can be partially specified using the three- 
jet variables described above. This is accomplished by first reducing the multijet system 
to a three-body system. A four-jet system is reduced to a three-body system by 
combining the two-jets with the lowest two-jet mass. The resulting three-body system 
can then be described using the variables 777,4 j, -'^s', -^^4', cos ^3/, ipy, fs'jf^', and f^i, 
where the primes remind us that two jets have been combined. A five-jet system is 
reduced to a three-body system by first combining the two-jets with the lowest two-jet 
mass to obtain a four-body system, and then combining the two bodies with the lowest 
two-body mass to obtain a three-body system. The resulting three-body system can 
then be described using the variables m^j, X3// , X4// , cos ^3// , ipy, /3",/4", and f^n. 

To complete the description of four-jet (five-jet) events we must now specify a 
further four (eight) variables that describe how the multijet system has been reduced 
to a three-body system. Consider first the step in which a four-body system is reduced 
to a three-body system. We label the two objects being combined A and B with 
Ea > Eb, where Ea and Eb are energies in the four-body rest- frame. To describe the 
(AB)-system we use the following four variables: 

(a) The normalized masses /a and /b: 



fA ^ 



rriA 



and /b = 



TUB 



(7) 



rriNj 



rriNj 
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where mp^j is the mass of the multijet system. 

(b) The two-body energy sharing variable X^, defined in the muhijet rest-frame as 
the fraction of the energy of the (AB)-system taken by A: 



Ea + Eb ^ ^ 



(c) The two-body angular variable V'ab' defined in the multijet rest-frame as the 
angle between (i) the plane containing the (y4i?)-system and the average beam 
direction, and (ii) the plane containing A and B. The prime reminds us that in 
order to define ijj'^^ we have combined two bodies to obtain the (AB)-system. 
Note that: 



_ (Pa X 'Pb) ■ (Pab X 'Pay) 

^osi^AB = ^-^^ (9) 

l-t^A X -Tb 1 1 r'AB X r'AV I 



For five-jet events we must also specify the step in which the five-jet system is reduced 
to a four-body system. We label the two jets that are combined C and D with Eq > E^, 
and use the four variables fc, fn, Xc, and V'ci?- 

In summary, a four-jet system is described using 12 variables: m^j, Xy, X4/, cos 6*3/, 
V's'! fs'y /4') /s'! Ia, Ib, Xa, and iP'ab- ^ five-jet system is described using 16 variables: 

'^SJ, -^3", -^4", COs6'3", ^/'3", /3", /4// , /s", /a', /b' , X^', i'A' b' ^ fc fo, Xq, and Ipcf). 

Note that following the convention of Ref. [1] the primes indicate which parameters 
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are defined after one or two steps in wliicli two objects liave been combined. 

5 Results 

Tfie (4N-4) multijet variables described in tfie previous sections provide a set of in- 
dependent parameters that span the multijet parameter space in the multijet rest 
frame. In the following the (4N-4) multijet distributions are compared with QCD 
and phase-space model predictions. All distributions are inclusive. If there are more 
than N jets in an event, the N highest Ex jets are used to define the multijet system. 
It should be noted that at fixed multijet mass the J^Et > 420 GeV, AR > 0.9, and 
Et > 20 GeV requirements place restrictions on the available multijet parameter space. 
Consequently, some regions of parameter space are depopulated due to a low exper- 
imental acceptance. These inefficient regions can be largely avoided in the three-jet 
analysis by placing suitable requirements on the multijet mass, leading-jet angle, and 
leading-jet Dalitz variable. In the following we have required that nisj > 600 GeV/c^, 
I cos ^3 1 < 0.6, and Xs < 0.9. These requirements select 1021 events with three-or- 
more jets, of which 320 events have more than three jets. Events entering the inclusive 
four-jet distributions are required to have m4j > 650 GeV/c^, |cos^3'| < 0.8, and 
X2' < 0.9. These requirements select 1273 events with four-or-more jets, of which 245 
events have more than four jets. Only 226 events enter into both the inclusive three- 
jet and inclusive four-jet distributions.. Note that the four-jet requirements reduce, 
but do not completely eliminate, the regions of low experimental acceptance. A more 
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restrictive X3/ requirement could be used to remove events populating the remaining 
region of low acceptance, but would cost a large reduction in statistics. Given the 
limited statistics of the present data sample, we have chosen to tolerate some regions 
of low experimental acceptance and use the phase-space model predictions to under- 
stand which regions of parameter space are affected. Finally, the inclusive five-jet 
data sample has very limited statistics, and we have therefore chosen to apply only 
the requirement m^j > 750 GeV/c^ to events entering the five-jet distributions. This 
requirement selects 817 events with five-or-more jets, of which 146 events have more 
than five jets. Only 148 events enter into both the five-jet and four-jet distributions, 
and only 42 events enter into both the five-jet and three-jet distributions. 

5.1 Multijet Mass Distributions 

In Ref. [0] HERWIG and NJETS QCD calculations were shown to give a good descrip- 
tion of the shapes of the observed multijet mass distributions for exclusive samples 
of high-mass multijet events. In Figs. ||a, |I]b, and |I]c both the HERWIG and NJETS 
predictions are shown to give good descriptions of the shapes of the inclusive m^j, m4j, 
and msj distributions for the high-mass multijet samples described in this paper. Note 
that over the limited mass range of the present data sample, to a good approximation 
the mjvj distributions are falling exponentially with increasing mass. 
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5.2 Three-Body Dalitz Distributions 

We begin by considering the inclusive three-jet Dahtz distributions. Event populations 
in the (X3, X4)-plane are shown in Fig. || for (a) data, (b) NJETS, (c) HERWIG, and 
(d) phase-space model predictions. The phase-space population is uniform over the 
kinematically allowed region. Neither the data nor the QCD predictions exhibit large 
density variations in the (X3, X4)-plane in the region of interest (X3 < 0.9), although 
with the relatively high statistical precision of the NJETS predictions the tendency for 
the predicted event density to increase as X4 becomes large is visible (note that as X4 
1 the third-to-leading jet Dalitz variable X5 — > 0). The observed X3 distribution is 
compared with phase-space model and QCD predictions in Fig. |a. The corresponding 
comparisons for the X4 distribution are shown in Fig. The HERWIG and NJETS 
predictions give reasonable descriptions of the observed distributions. Note that the 
observed distributions are not very different from the phase-space model predictions. 

We next consider the inclusive four-jet distributions. Event populations in the 
(X3/, X4/)-plane are shown in Fig. | for (a) data, (b) NJETS, (c) HERWIG, and 
(d) phase-space model predictions. The phase-space population is not uniform over 
the kinematically allowed region. Care must therefore be taken in interpreting the 
distributions. The data and the QCD predictions exhibit a more uniform event density 
over the {Xy, X4/)-plane. The observed X3/ distribution is compared with phase- 
space model and QCD predictions in Fig. The corresponding comparisons for 
the X4/ distribution are shown in Fig. |^. The HERWIG and NJETS predictions 

19 



give reasonable descriptions of the observed distributions. Note that compared to the 
phase-space model predictions, the data and QCD predictions prefer topologies with 
large X31 and large X4/ (note that as Xy 1 the three-body topology approaches a 
two-body configuration, and as X41 — >• 1 we have X^i —>■ 0). 

Finally, consider the inclusive five-jet distributions. Event populations in the (X3//, 
X4//)-plane are shown in Fig. ^ for (a) data, (b) NJETS, (c) HERWIG, and (d) phase- 
space model predictions. Again, the phase-space population is not uniform over the 
kinematically allowed region, and care must be taken in interpreting the distributions. 
The observed event population and the QCD predictions are more uniformly distributed 
over the (X3//, X4//)-plane. However, all distributions are depleted as X-y — > 1 and 
X4// 1. The observed Xy distribution is compared with phase-space model and QCD 
predictions in Fig. |^. The corresponding comparisons for the X41/ distribution are 
shown in Fig. |^. The HERWIG and NJETS predictions give reasonable descriptions of 
the observed distributions. Note that compared to the phase-space model predictions, 
the data and QCD predictions prefer topologies with X-y 1 and X4// — > 1. 

5.3 Three-Body Angular Distributions 

We begin by considering the inclusive three-jet angular distributions. Event popula- 
tions in the (cos6'3, ip3)-p\a.ne are shown in Fig. ^ for (a) data, (b) NJETS, (c) HERWIG, 
and (d) phase-space model predictions. The phase-space population is approximately 
uniform. In contrast both the observed distribution and the QCD predictions exhibit 
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large density variations over the (cos6'3, ■?/'3)-plane, with the event density increasing 
as I cos^^sl 1 and ip^ ^ or tt . The increase in event rate as | cos^^al ^ 1 is similar 
to the behaviour of the leading-jet angular distribution resulting from the 2 2 LO 
QCD matrix element. The increase in event rate as t/'s — > or tt reflects the prefer- 
ence of the three-jet matrix element for topologies which are planar. It is interesting 
to note that as cos^^s — > 1 the NJETS calculation shows a preference for configura- 
tions with ip^ ^ rather than vr and as cos^^3 ^ — 1 the NJETS calculation shows a 
preference for configurations with ip^ ^ n rather than 0. These preferred regions of 
the parameter space correspond to configurations in which jet 5 is closer to the beam 
direction, and therefore reflect the initial state radiation pole in the matrix element. 
The observed cos 63 distribution is compared with phase-space model and QCD predic- 
tions in Fig. I^a. The corresponding comparisons for the ^3 distribution are shown in 
Fig. Both HERWIG and NJETS predictions give reasonable descriptions of the ob- 
served distributions, which are very different from the phase-space model predictions. 
Note that the observed cos ^3 distribution is also very similar to the LO prediction for 
qq qq scattering [l^ . 



Next, consider the inclusive four-jet angular distributions. Event populations in 



the (cos ^3/, ^3/)-plane are shown in Fig. |10| for (a) data, (b) NJETS, (c) HERWIG, 



and (d) phase-space model predictions. The phase-space population is approximately 
uniform. In contrast both the observed distribution and the QCD predictions exhibit 
large density variations over the (cos6'3/, ^/'3/)-plane, with the event density increasing 
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as I cos^^s'l — > 1 and ipy ^ or tt. This behaviour is similar to the behaviour of the cor- 
responding three-jet distributions. The observed cos 6*3/ distribution is compared with 



phase-space model and QCD predictions in Fig. 11a. The corresponding comparisons 



for the ip^i distribution are shown in Fig. |TT]b. Both HERWIG and NJETS predictions 



give reasonable descriptions of the observed distributions, which are very different from 
the phase-space model predictions. Note that the observed cos ^3/ distribution is also 
very similar to the LO prediction for qq qq scattering. 

Finally, consider the inclusive five-jet angular distributions. Event populations in 



the (cos%., ^/'g'O-plane are shown in Fig. |T| for (a) data, (b) NJETS, (c) HERWIG, 



and (d) phase-space model predictions. The phase-space population is not uniform, 
and care must therefore be taken in interpreting the distributions. However, both the 
observed distribution and the QCD predictions exhibit much larger density variations 
over the (cos^/, '?/'3//)-plane, with the event density increasing as |cos03"| 1 and 
ip^j, — >^ or TT. This behaviour is similar to the behaviour of the corresponding three-jet 
distributions. The observed cos ^3" distribution is compared with phase-space model 



and QCD predictions in Fig. 13a. The corresponding comparisons for the ifj^," distri- 



bution are shown in Fig. |13|b. Both HERWIG and NJETS predictions give reasonable 
descriptions of the observed distributions, which are very different from the phase-space 
model predictions. Note that the observed cos 6*3/' distribution is also very similar to 
the LO prediction for qq —>■ qq scattering even though there are now five jets in the 
final state. 
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5.4 Single-Body Mass Distributions for Three-Body Systems 



The single-jet mass fraction distributions are shown in Fig. [1^ for inclusive three-jet 
events. The fj distributions are reasonably well described by the HERWIG Monte 
Carlo predictions, although there is a tendency for the HERWIG fragmentation model 
to slightly overestimate the fraction of low-mass jets. The observed distributions peak 
at fj ~ 0.05 or less. Hence, for many purposes, jets at high energy can be considered 
to be massless. Note that since jets are massless in the matrix element calculations, 
there are no NJETS predictions for the fj distributions. 

The /-,' and fj" distributions are shown for inclusive four-jet and inclusive five-jet 



events in Figs. |T5| and ^ respectively. These distributions exhibit a single-jet peak 
at low mass-fractions (less than 0.05), and have a long tail associated with two-jet j' 
systems, and two-jet or three-jet j" systems. The HERWIG predictions give a good 
description of all the distributions except perhaps at very low mass fractions (less than 
0.05) where there is tendency to overestimate the observed jet rate. Although the 
NJETS calculations do not provide predictions for the single-jet part of the fj' and fj" 
distributions, they are seen to correctly predict the tail associated with multijet j' and 
j" systems. 

5.5 Two-Body Energy Sharing Distributions 

The observed Xa distribution is shown in Figs. |l7|a and to be reasonably well 
described by the HERWIG and NJETS predictions. The data and the QCD predictions 
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favor a more asymmetric sharing of energy between the two jets A and B than predicted 
by the phase-space modeL This reflects the presence of the soft gluon radiation pole 
in the QCD matrix element. In Figs. p!7|c and Od the Xa' distributions are shown 



to be quahtatively similar to the corresponding Xa distributions, and also similar to 
the corresponding Xc distributions shown in Figs. p!7| e and ^Llf . In general the data 



are reasonably well described by the QCD predictions and are very different from the 
phase-space model predictions. 

5.6 Two-Body Angular Distributions 

The observed iP'ab distribution is shown in Figs. ISa and |I3b to be well described 



by the HERWIG and NJETS predictions. The phase-space model prediction is also 
approximately uniform, but underestimates the fraction of events in which the plane 
of the two-body system is close to the plane containing the two-body system and the 
beam direction {^iPab ^ or tt). In Figs. |l8|c and pTSjd the V'a'B' distributions are shown 



to be qualitively similar to the corresponding ip'j^^ distributions. The ipcn distributions 



shown in Figs. 18 e and 18 f are similar to the phase-space model predictions. In all cases 
the data are well described by the QCD predictions. None of the observed distributions 
are very different from the phase-space model predictions, although the phase-space 
model calculation does underestimate the event rate as iP'ab — > or tt, or as iPa'b' 

or TT. 
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5.7 Single-Body Mass Distributions for Two-Body Systems 

The observed /a, fs, /a', Ib', fc, and distributions are shown in Figs. |T^, |T^, 
IP|c, P^d, and respectively to be reasonably well described by the HERWIG 
predictions although there is a tendency for the HERWIG predictions to overestimate 
the jet rate at very small single-jet masses. In all cases the distributions exhibit a 
single-jet mass peak at small mass fractions (~ 0.02 or less). The Ja' and Jb' have a 
long high-mass tail which corresponds to two-jet A' and B' systems. This tail is well 
described by the NJETS predictions. 

5.8 Test 

In general both NJETS and HERWIG predictions give a good first description of the 
observed multijet distributions, which correspond to (4N-4) variables that span the 
N-body parameter space. A more quantitative assessment can be made by examining 
the psr degree of freedom that characterizes the agreement between the observed 
distributions and the QCD predictions. The are listed for each distribution in Table 
1. The computed x^'s take into account statistical uncertainties on both measured 
points and the QCD Monte Carlo predictions, but do not take into account systematic 
uncertainties. In Ref. ^ the systematic uncertainties were mapped out for a limited 
set of multijet distributions, and found to be small compared to statistical uncertain- 
ties. Unfortunately, for the more complicated multijet parameter space of the present 
analysis, limited computing resources do not allow us to fully map out the systematic 
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uncertainties on the predictions. However, even in the absence of a full evaluation of 
the systematic uncertainties, an examination of Table 1 shows that NJETS provides a 
reasonable description of all of the observed multijet distributions except perhaps the 
Xa distribution. The combined for the NJETS description of all of the three-jet 
distributions x^/NDF = 1.03 (45 degrees of freedom). The corresponding result for 
the four-jet distributions is x^/NDF = 1.47 (63 degrees of freedom) if the Xa distri- 
bution is included in the comparison, and x^/NDF = 1.24 (55 degrees of freedom) if 
the Xa distribution is not included. The result for the combined five-jet distributions 
is x^/NDF = 1.21 (63 degrees of freedom). The observed distributions are described 
less well by the HERWIG parton shower Monte Carlo predictions, for which the X4, 
cos 6*3/, ■1/^3/, and cos 6*3" distributions have x^s significantly poorer than those for the 
corresponding NJETS predictions. Restricting the comparison to those distributions 
predicted by both the NJETS and HERWIG calculations (i.e. all distributions except 
the single-body mass fraction distributions) we find the overall x^ per degree of free- 
dom for the HERWIG comparison of the combined three-jet distributions is x^/NDF 
= 1.58 (45 degrees of freedom), for the combined four-jet distributions x^/NDF = 1.63 
(63 degrees of freedom), and for the combined five-jet distributions x^/NDF = 1.52 
(63 degrees of freedom). 

6 Conclusions 

The properties of high-mass three-jet, four-jet, and five-jet events produced at the Fer- 
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milab Tevatron proton-antiproton collider have been compared with NJETS LO QCD 
matrix element predictions, HERWIG QCD parton shower Monte Carlo predictions, 
and predictions from a model in which events are distributed uniformly over the avail- 
able multibody phase-space. The phase-space model is unable to describe the shapes 
of multijet distributions in regions of parameter space where the QCD calculations pre- 
dict large contributions from initial- and final-state gluon radiation. In contrast, the 
QCD predictions give a good first description of the observed multijet distributions, 
which correspond to (4N-4) variables that span the N-body parameter space. A more 
quantitative assessment based on the per degree of freedom that characterizes the 
agreement between the observed distributions and the QCD predictions shows that 
NJETS gives a good description of all the distributions except perhaps the Xa distri- 
bution for four-jet events. The NJETS predictions seem to give a better description 
of the observed distributions than the HERWIG predictions. This is particularly true 
of the X4, COSTS', -03', and cos ^3" distributions. Finally, we do not see clear evidence 
for any deviation from the predicted multijet distributions that might indicate new 
phenomena associated with the presence of many hard partons in the final state. The 
general agreement between data, NJETS, and HERWIG suggests that 2 — > 2 scattering 
plus gluon radiation provides a good first approximation to the full LO QCD matrix 
element for events with three, four, or even five jets in the final state. 
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Table 1: Statistical comparison of agreement between oberserved and predicted dis- 
tributions. The per degree of freedom are listed for comparisons of the various 
observed and QCD predicted distributions shown in the figures. 
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Figure 1: Inclusive multijet mass distributions for topologies with (a) three jets, (b) 
four jets, and (c) five jets. Observed distributions (points) are compared with HERWIG 
predictions (triangles) and NJETS predictions (squares). 
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Figure 2: Three-jet Dalitz distributions after imposing the requirements msj > 
600 GeV/c2, Xs < 0.9, and | cosOs] < 0.6, shown for (a) data, (b) NJETS, (c) HER- 
WIG, and (d) the phase-space model. 
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Figure 3: Inclusive three-jet Dalitz distributions for events that satisfy the requirements 
TTisj > 600 GeV/c^, X3 < 0.9, and | costal < 0.6. Data (points) are compared with 
HERWIG predictions (triangles), NJETS predictions (squares), and phase-space model 
predictions (curves) for (a) X^, and (b) X4. 
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Figure 4: Inclusive four-jet Dalitz distributions for events that that satisfy the require- 
ments m4j > 650 GeV/c^, X^i < 0.9, and |cos^3/| < 0.8, shown for (a) data, (b) 
NJETS, (c) HERWIG, and (d) phase-space model predictions. 
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Figure 5: Dalitz distributions for inclusive four-jet topologies that satisfy the require- 
ments m4j > 650 GeV/c^, X3/ < 0.9, and | cos6'3/| < 0.8. Data (points) are compared 
with HERWIG predictions (triangles), NJETS predictions (squares), and phase-space 
model predictions (curves) for (a) X^i, and (b) X4/. 
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Figure 6: Inclusive five-jet Dalitz distributions for events that satisfy the requirement 
msj > 750 GeV/c2, shown for (a) data, (b) NJETS, (c) HERWIG, and (d) phase-space 
model predictions. 
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Figure 7: Dalitz distributions for inclusive five-jet topologies that satisfy the require- 
ment msj > 750 GeV/c^. Data (points) are compared with HERWIG predictions 
(triangles), NJETS predictions (squares), and phase-space model predictions (curves) 
for (a) X3//, and (b) X4//. 
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Figure 8: Inclusive three-jet angular distributions for events that satisfy the require- 
ments msj > 600 GeV/c^, X3 < 0.9, and |cos6'3| < 0.6. Event populations in the 
(008 6*3, ip3)-plane are shown for (a) data, (b) NJETS, (c) HERWIG, and (d) phase- 
space model predictions. 
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Figure 9: Inclusive three-jet angular distributions for events that satisfy the require- 
ments msj > 600 GeV/c^, X3 < 0.9, and | cos 6*3! < 0.6. Data (points) are compared 
with HERWIG predictions (triangles), NJETS predictions (squares), and phase-space 
model predictions (curves) for (a) cos ^3 and (b) ■03- The broken curve in the cos ^3 
figure is the LO QCD prediction for qq — > qq scattering. 
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Figure 10: Inclusive four-jet angular distributions for events that satisfy the require- 
ments m4j > 650 GeV/c^, X^/ < 0.9, and |cos6'3'| < 0.8. Event populations in the 
(008 6*3', ■V'3')-plane are shown for (a) data, (b) NJETS, (c) HERWIG, and (d) phase- 
space model predictions. 
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Figure 11: Inclusive four-jet angular distributions for events that satisfy the require- 
ments m^j > 650 GeV/c^, X3/ < 0.9, and | cos 6*3/! < 0.8. Data (points) are compared 
with HERWIG predictions (triangles), NJETS predictions (squares), and phase-space 
model predictions (curves) for (a) costs' and (b) ■03'- The broken curve in the costs' 
figure is the LO QCD prediction for qq — > qq scattering. 
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Figure 12: Inclusive five-jet angular distributions for events that satisfy the requirement 

m^j > 750 GeV/c^. Event populations in the (cos ^3", V's")"?!^'^^ are shown for (a) 
data, (b) NJETS, (c) HERWIG, and (d) phase-space model predictions. 
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Figure 13: Inclusive five-jet angular distributions for events that satisfy the requirement 
TTT'^j > 750 GeV/c^. Data (points) are compared with HERWIG predictions (triangles), 
NJETS predictions (squares), and phase-space model predictions (curves) for (a) cos 
and (b) ■03"- The broken curve in the costs'/ figure is the LO QCD prediction for 
qq — > qq scattering. 
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Figure 14: Single-jet mass fraction distributions for inclusive three-jet events. Data 
(points) compared with HERWIG predictions (triangles), shown for (a) the highest 
energy jet in the three-jet rest-frame, (b) the second-to-highest energy jet, and (c) the 
third-to-highest energy jet. 
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Figure 15: Single-body mass fraction distributions for inclusive four-jet events. Data 
(points) compared with HERWIG predictions (triangles), and NJETS predictions (his- 
tograms), shown for (a) the highest energy body in the three-body rest-frame, (b) the 
second-to-highest energy body, and (c) the third-to-highest energy body. 
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Figure 16: Single-body mass fraction distributions for inclusive five-jet events. Data 
(points) compared with HERWIG predictions (triangles), and NJETS predictions (his- 
tograms), shown for (a) the highest energy body in the three-body rest-frame, (b) the 
second-to-highest energy body, and (c) the third-to-highest energy body. 
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Figure 17: The two-body energy sharing distributions for inclusive four-jet and five-jet 
events. Data (points) are compared with HERWIG predictions (triangles), NJETS 
predictions (squares), and phase-space predictions (curves) for (a) Xa, (b) after 
dividing by the phase-space model predictions, (c) Xa', (d) Xa' after dividing by the 
phase-space model predictions, (e) Xc, and (f) Xc after dividing by the phase-space 
model predictions. 
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Figure 18: Two-body angular distributions for inclusive four-jet and five-jet events. 
Data (points) are compared with HERWIG predictions (triangles) , N JETS predictions 
(squares), and phase-space predictions (curves) for (a) V'ab; (b) i^AB after dividing by 
the phase-space model predictions, (c) iI>%b'-i (d) ^a'B' after dividing by the phase- 
space model predictions, (e) V'c'D) and (f) -0^^ after dividing by the phase-space model 
predictions. 
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Figure 19: Single-body mass fraction distributions for two-body systems in inclusive 
four-jet and five-jet events. Data (points) are compared with HERWIG predictions 
(triangles), and NJETS predictions (histograms) for (a) /a, (b) /b, (c) /a', (d) /g/, (e) 
fc, and (f) fo- 
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